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SUMMARY 


A theoretical method has been developed for computing approximate laminar 
heating rates on three-dimensional configurations at angle of attack. The 
method is based on the axisymmetric analogue which is used to reduce the three- 
dimensional boundary-layer equations along surface streamlines to an equivalent 
axisymmetric form by using the metric coefficient which describes streamline 
divergence (or convergence) . The method has been coupled with a three- 
dimensional inviscid flow field program for computing surface streamline paths, 
metric coefficients, and boundary-layer edge conditions. Good agreement has 
been obtained with experimental data on a spherically blunted 15° half-angle 
cone, a spherically blunted 80° sweep slab delta wing, and a Space Shuttle 
Orbiter type configuration at angles of attack up to 25°. The method provides 
a useful means of computing heating rates on advanced entry configurations. 


INTRODUCTION 

With the development of the Space Shuttle transportation system well under 
way, attention is being focused on more advanced Earth orbital transportation 
systems which will be expected to replace the Space Shuttle system near the end 
of this century. These vehicles will be subjected to severe aerodynamic heating 
during reentry, and the design of their thermal protection systems will require 
detailed information on surface heating rates. Three-dimensional winged config- 
urations are currently being considered (refs. 1 and 2) and they will probably 
fly at relatively large angles of attack. It will first be necessary to compute 
the three-dimensional viscous flow over the vehicle in order to accurately pre- 
dict the surface heating rates. 

With advances in high-speed computers and numerical techniques, much prog- 
ress has been made in recent years in computing viscous compressible flows using 
the Navier-Stokes equations (ref. 3). However, solutions for general three- 
dimensional winged configurations are far beyond the "state of the art" of cur- 
rent applications, and their demand on computer resources would far exceed the 
capabilities of the current generation of computers. 

Much of the needed information on surface heating rates can be obtained 
from a "classical approach" where the outer inviscid flow field is computed 
independent of the boundary layer and is used to provide edge conditions for a 
three-dimensional boundary-layer calculation near the surface. Solutions have 
been obtained for the general three-dimensional compressible boundary-layer 
equations. (See, for example, refs. 4 to 8.) To date, the application of these 
methods has been limited to relatively simple geometries. An attempt has been 
made to extend the method of reference 5 to more realistic configurations; but, 
it has not yet been successful. However, even if these methods could be applied 
to the configurations of interest in the present study, their use would be cost 
prohibitive for design calculations. 



Cooke has developed an axi symmetric analogue for the general three- 
dimensional compressible boundary- layer equations (ref. 9) which greatly simpli- 
fies the calculation of three-dimensional viscous flows. Following that 
approach, the general three-dimensional boundary- layer equations are written in 
streamline coordinates and the cross-flow velocities (tangent to the surface and 
normal to streamline direction) are assumed zero. This reduces the three- 
dimensional boundary-layer equations to a form identical to that for axisymmet- 
ric flow provided: (1) the distance along a streamline is interpreted as dis- 

tance along an "equivalent axi symmetric body" and (2) the metric that describes 
the spreading of the streamlines is interpreted as the radius of the equivalent 
axisymmetric body. This allows any existing axisymmetric boundary-layer program 
to be used to compute three-dimensional heating rates along a streamline in 
regions where the small cross-flow assumption is valid. By considering multiple 
streamline paths, an entire vehicle can be covered. 

Hayes (ref. 10) has shown that the cross flow in the boundary layer is 
small when the streamline curvature is small. Vaglio-Laurin (ref. 11) has shown 
that when the wall is highly cooled the cross flow in the boundary layer is 
small even for cases where the streamline curvature is not small. Further, in 
references 12 to 14, for example, comparisons with experimental data and other 
theoretical calculations, including some cases where the cross flow in the 
boundary layer is not small, indicate that reasonably accurate heating rates can 
be obtained using the axisymmetric analogue. 

The most obvious difficulty in applying this technique is that of computing 
the three-dimensional inviscid flow field from which the surface streamline pat- 
tern and the metrics are determined. In reference 15, DeJarnette and Hamilton 
present a relatively simple method of calculating the streamline location and 
metric coefficient along the streamline from a known surface pressure distribu- 
tion. With this approach the path of a single streamline is computed on the 
surface; simultaneously, a second-order differential equation is solved along 
the streamline for the metric. The major difficulty with this approach is that 
second derivatives of the surface pressure are required to calculate the metric. 
These are very difficult to obtain with sufficient accuracy from either experi- 
mental measurements or approximate methods (such as Newtonian theory) except for 
a few simple cases. In reference 16, DeJarnette presents an improved method for 
calculating the inviscid surface streamlines. Two adjacent streamlines are com- 
puted from a known surface pressure distribution, using only first derivatives 
of the pressure; and the metric is obtained by determining the rate at which the 
streamlines spread on the surface. This is a much better approach since it 
relaxes some requirements upon the accuracy of the input surface pressure dis- 
tribution. However, it still requires more detailed pressure distributions than 
are generally available from most experimental studies. 

In both references 17 and 18, methods have been presented for computing the 
metric coefficient on the surface of a body from an inviscid flow field solu- 
tion. However, in each case it would be very difficult to accurately compute 
the metric coefficient at the boundary-layer edge using these methods. Although 
this will not be considered in the present paper, it is important that the cho- 
sen technique should not have such a limitation. Since in the present approach, 
the metric coefficient is computed from information on two adjacent streamlines. 


2 


■■I I HIM II ■■■MU 



it will be very easy to extend the method to use boundary-layer edge conditions 
to compute the metric coefficient at the boundary-layer edge. 

Recently, a method has been developed for computing the three-dimensional 
inviscid flow field over general configurations at angle of attack. (See 
refs. 17 and 19.) Using information about the velocity directions from the 
inviscid solution, streamlines can be traced on the vehicle surface. Then, a 
method similar to that presented by DeJarnette (ref. 16) can be used to compute 
the metrics along each streamline. This approach offers the advantage of com- 
puting heating rates on three-dimensional vehicles without requiring experimen- 
tal pressure distributions. 

In the present paper, an axisymmetric analogue method for computing heating 
rates on three-dimensional configurations at angle of attack is developed. The 
necessary inviscid streamline information was obtained using the inviscid- 
solution procedure of reference 17. The accuracy of this solution technique is 
demonstrated by direct comparison of computed heating rates with experimental 
data. 


The information presented in this report was included in a report entitled 
"Calculation of Heating Rates on Three-Dimensional Configurations" submitted in 
partial fulfillment of the requirements for the Degree of Engineer, The George 
Washington University, December 1979. 

SYMBOLS 

A parameter in transformation from cylindrical to local polar 

coordinates 

A parameter in transformation from local polar to cylindrical 

coordinates 

Cp specific heat at constant pressure, J/kg-K 

c v specific heat at constant volume, J/kg-K 

h static enthalpy, J/kg 

h e local enthalpy at edge of boundary layer, J/kg 

h r reference heat-transf er coefficient, W/m 2 -s-K 

h s stagnation-point heat-transfer coefficient, W/m 2 -s-K 

h w heat-transfer coefficient at wall, W/m 2 -s-K 

h-j metric coefficient in £ -direction, m 

h 2 metric coefficient in g-direction, m 
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total enthalpy, J/kg 
thermal conductivity, W/m-K 
streamline curvature, 1/m 
length, m 

distance measured along leading edge of delta wing, m 
Mach number 

Prandtl number, Cpjj/k 
unit Reynolds number 
distance normal to surface, m 
pressure, N/m^ 
heat-transfer rate, W/m^ 

distance normal to streamline on surface, m 
recovery factor 

radius of curvature of spherical nose, m 

displacement of coordinate pole, m 

local polar coordinates (see fig. 4) 

cylindrical coordinates (see fig. 4) 

constant in Sutherland's viscosity law, K 

distance along streamline, m 

arc length, m 

temperature, K 

ratio of enthalpies, h e /H e 

velocity component in ^-direction, n\/s 

free-stream velocity, m/s 

velocity component in 3-direction, m/s 

velocity component defined by equation (B1 ) , m/s 
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v n velocity component defined by equation (B2) , m/s 

w velocity component in n-direction, m/s 

w t velocity component defined by equation (B3) , m/s 

w-j f W 2 /W 3 Cartesian velocity components defined by sketch (a) of 
appendix B, m/s 

x,y, z Cartesian coordinates (see fig. 5) 

a angle of attack 

3 coordinate normal to streamline on surface 

3 boundary-layer velocity gradient parameter 

T geometry angle defined by equation (22) (see sketch (c) of appendix B) 

Y ratio of specific heats, c p /c v 

6^ geometry angle defined by equation (21) (see sketch (b) of appendix B) 

C w ratio of total enthalpies, H w /H e 

C ^ slope of total enthalpy profile at wall, 

H transformed boundary- layer coordinate normal to wall 

0 angle defining velocity direction on surface measured relative to 

meridional plane (see fig. 8) 

M viscosity, N-s/m 2 

£ coordinate along streamline on surface 

Z transformed boundary- layer coordinate along surface 

streamline 

p density, kg/m^ 

$ dissipation function, W/m 3 

Subscripts: 
b body 

max maximum 
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e edge of boundary layer 

s inviscid stagnation point 

w wall 

00 free stream 


METHOD OF SOLUTION 

In the sections that follow, a method for calculating heating rates on 
three-dimensional configurations using the axisymmetric analogue (ref. 9) is 
described, and the equations necessary for applying the technique are developed. 
The assumptions made from the outset are that: (1) the Reynolds number of the 

flow is sufficiently large so that boundary-layer concepts apply; (2) the fluid 
is a perfect gas; (3) the flow is laminar; and (4) the entropy at the edge of 
the boundary layer is constant. Although the last three assumptions are not 
necessary, they will simplify the development of the technique. 


Axisymmetric Analogue 

Consider an orthogonal streamline coordinate system (£,g,n) where £ is 
directed along an inviscid surface streamline and n is distance measured along 
the outward normal from the surface (fig. 1). The element of arc length di in 
such a curvilinear coordinate system is given by the equation 


ds 2 


= h-| d£ 2 + h 2 <33 2 + dn 2 


( 1 ) 


where h] and h 2 are metric coefficients. By restricting attention to a thin 
boundary-layer region near the surface, both h^ and h 2 can be assumed to be 
functions of £ and 8 only. 

Taking the velocity components in the 8-, and n-directions as u, v, 
and w, respectively, the general three-dimensional boundary-layer equations can 
be written as follows (ref. 9): 


Continuity 


1 

h 1 h 2 


9 

— (Ph 2 u) 


9 

+ — - (Ph-j v) 


9 

+ — (ph-|h 2 w) 
9n 


0 


( 2 ) 
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g -momentum 


( u 3u v 9u 3u uv 

+ + w — + 

hi 3£ h 2 33 3n ^1^2 ^3 


1 


v 2 Sh^ 
h-|h 2 3£ 


1 3p 3 / 3u' 

h-| 3£ 3n\ 3n/ 


3 -momentum 


/u 3v v 3v 3v uv ^2 

lh] 3£ h 2 36 3n h]h 2 3£ 


u 


2 3hi 


h 1 h 2 33 


1 3p 3 / 3v^ 

h 2 33 3n\ 3ny 


n-momentum 


3p 

3n 


= 0 


(3) 


(4) 


(5) 


Energy 


/ u 3T v 

P c p — — + — 
p lhi 3£ h 2 


3t 3t\ 

— + w — = $ 

33 3n ) 


u 3p 


+ h-j 3£ 


v 3p 3p 3 / 3 t\ 

+ — — + w — + — k — 

h 2 33 3n 3n y 3n J 


( 6 ) 


where $ is the dissipation function. (See ref. 9.) 

It is convenient to express derivatives in the ^-direction in terms of the 
distance s along a streamline as 

1 3 _ 3 

h] 3 Z, 3 s 
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Further, it is assumed that the cross-flow velocity v and its derivatives are 
small. With these simplifications, equations (2), (3), (5), (6), and (4) are as 
follows : 


Continuity 


3 3 

— (Ph 2 u) + — (ph 2 w) = 0 
3s 3n 


s-momentum 


/ 3u 

3u\ 

3 P 

3 / 

3u 

P u — 

+ w — = 

— 

+ — n 

— 

\ 3s 

3n/ 

3s 

3n\ 

3n 


n-momentum 


3p 

3n 


0 


Energy 


Pc p 


3t 

— + w 

3s 



9p 3 / 3 t 

u — + — k — 

3s 3n\ 3n 


/3u 

\3n 


3 -momentum 


/ 3v 3v uv 9 ^2 

p u — + w — + — 

l 3s 3n h 2 3s 



1 3p 3 / 3v\ 

h 2 33 3n\ 3n/ 


where k is the streamline curvature term given by the equation 


1 9h l 

~ h-|h 2 33 


(7) 


( 8 ) 


(9) 


( 10 ) 


(ID 


( 12 ) 
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Equations (7) to (10) are identical to the usual axisymmetric boundary- layer 
equations (ref, 20) if s is interpreted as distance along an "equivalent 
axisymmetric body" and h 2 is interpreted as the radius of that body. This 
set of equations can be solved for u, w, P, and T through the boundary 
layer along a streamline using one of the many methods that have been developed 
for solving the axisymmetric boundary-layer equations (for example, refs. 21 
and 22) . With u, w, p, and T assumed known from a previous solution, equa- 
tion (11) becomes a linear equation which can easily be solved for v, the 
cross-flow velocity component. In the present application, only the heat trans- 
fer is of interest; thus it is unnecessary to solve the cross-flow momentum 
equation since in the axisymmetric analogue it has no effect on heating rates. 


Coordinate System and Surface Geometry 

The present paper is concerned with calculating heating rates on three- 
dimensional configurations - typical of the configuration illustrated in fig- 
ure 2 - which are symmetric about the XZ-plane (pitch plane). The fuselage of 
these vehicles can be easily described in cylindrical coordinates (r b , 4>,x) 
using a function of the form (fig. 3) 


r b “ 


( 13 ) 


where the pole of the coordinate system is located on the X-axis. However, when 
wings extend from the fuselage, equation (13) may become multivalued. (See 
fig. 4.) One solution to this problem is to allow the pole to be a function 
of x 


cq = r G< x ) 


(14) 


Then, at a given x-station, the pole can be located such that r b is a single- 
valued function. This new coordinate system R b ,$,x will be called a "local 
polar coordinate system" in which the vehicle surface can be represented by a 
function of the form 


Rb = «b (*' X) 


( 15 ) 


which is much more general than equation (13). It should be noted that the new 
coordinate system is generally nonorthogonal; however, for the special case 
where rg(x) = 0, it reduces to the standard cylindrical coordinate system which 
is orthogonal. 

Vachris and Yaeger (ref. 23) have developed a computer program called QUICK 
which uses the local polar coordinate system to describe vehicle geometries. In 
this program, analytic curves are patched together to give a functional repre- 
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sentation of the vehicle geometry similar in form to equation (15). The QUICK 
program is compatible with the method that will be used to obtain the inviscid 
flow field solution (refs. 17 and 19). Thus, it will be used exclusively in the 
present study to describe vehicle surface geometry. 

The detailed relationship between local polar and cylindrical coordinates 
on a surface are given in appendix A, but the transformation equations are sum- 
marized here for convenience. First, the equations describing the transforma- 
tion from cylindrical coordinates rb,<i>,x to local polar coordinates Rfc,,$,X 
on a surface are 


Rb = r^A 2 + 1 cos <(> (16a) 

$ = tan - ^ A (16b) 


X = x 


0 6 c) 


where 


A = tan <}> - r 0 (x)/(rk cos 4>) (16d) 

Similarly, the inverse transformation from local polar coordinates Rb,$,X to 
cylindrical coordinates on a surface is expressed as 


r b = R b^A 2 + 1 cos $ 

(17a) 

<p = tan " 1 A 

(17b) 

x = X 

(17c) 


where 


A = tan $ + ro(X)/(Rb cos $) 


(1 7d) 


Since the vehicle geometry will be obtained from the QUICK program in the func- 
tional form of equation (15), the inverse transformation equations should be 
more useful. 
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Calculation of Inviscid Surface Streamlines 


and Metric Coefficients 

One of the most difficult steps in applying the axisymmetric analogue is 
calculating the inviscid surface streamlines and the metric coefficient h 2 
along each streamline. In the present study the inviscid streamline information 
will be calculated using the results of a set of three-dimensional inviscid 
flow field programs described in references 17, 19, and 24. The manner in which 
the flow field is calculated is illustrated in figure 5. First, the subsonic/ 
transonic flow over the nose is computed using a time dependent method BLUNT 
similar to that described in reference 24. Fran this solution an initial data 
plane (IDP) normal to the X-axis is obtained, in a region where the flow is 
entirely supersonic. Then, a finite-difference marching technique STEIN 
(refs. 17 and 19) is used to continue the solution downstream in the supersonic 
region step by step in planes normal to the X" axis * These solutions provide 
the inviscid flow field data over the entire vehicle which will be used to cal- 
culate the streamlines and metric coefficients. 

Nose region .- In the nose region of a blunt body, the surface streamlines 
originate at the stagnation point as illustrated in figure 6. Letting £ be 
the coordinate in the streamline direction and 3 the coordinate in a direction 
normal to the streamline and tangent to the surface, the differentials of arc 
length are 


ds = h] d£ (1 8a) 

and 

d q = h 2 d3 (1 8b) 


respectively. The metric coefficient h 2 gives a measure of the divergence or 
convergence of the streamlines on the surface and is the "equivalent radius" 
that must be used in the axisymmetric boundary-layer equations. 

The nose geanetry of most nonablating blunt reentry vehicles is smooth; 
thus, the pressure distribution in the nose region of these vehicles is also 
relatively smooth. Because of this, the pressures can be canputed by BLUNT and 
used as inputs to the computer program described in reference 1 6 to calculate 
the streamlines and metric coefficients fran the stagnation point to the station 
where the initial data plane intersects the body surface. (See fig. 6.) In 
principle, the same technique could be used to continue the streamline solutions 
downstream; but it could lead to large inaccuracies since the pressure distri- 
bution can change radically downstream of the nose on general three-dimensional 
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configurations. Thus, downstream of the IDP, a new technique for obtaining 
streamline solutions has been developed, and it will be described in detail in 
the next section. 

Downstream s u personic region .- DeJarnette (ref. 1 6) presents the following 
two equations which can be used to calculate the streamline position: 



•j /tan 0 cos 6^ 
rt> \ cos T 


- tan T sin 6 


and the metric coefficient h 2 

ft*) ./- litL 

Wx \r b cos 0/ 


( 19 ) 


( 20 ) 


The angles 6^ and T are the body-geometry angles defined in reference 15. 
(See also fig. 7.) They are given by the equations 


tan 



( 21 ) 


and 


tan T = cos 



( 22 ) 


The angle 0 is related to the velocity direction on the surface (fig. 8) and 
can be obtained from the inviscid flow field solution. 

The distance along a streamline can be calculated from the following dif- 
ferential equation: 


9s\ _ 1 

9x^13 cos 0 cos T 


(23) 
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Since 3 is constant on a given streamline, equation (19) can be inte- 
grated along a streamline to determine the variation of the cylindrical 
angle <(> with respect to x. Similarly, equation (23) can be integrated to 
give the streamline length. To obtain the metric coefficient along a stream- 
line, the primary streamline - the streamline along which heating rates are 
calculated - and an adjacent streamline (separated from the primary streamline 
by A3) are computed simultaneously. Then, at a given x-station the metric 
coefficient h 2 is computed from equation (20) after first approximating the 
derivative (9<t>/33) x numerically using one-sided differences. In reference 16 
numerical experiments were performed using several types of difference schemes 
to approximate (3<|>/93) x and the preceding procedure was found to be suffi- 
ciently accurate. 

The angles 6(j> and T given by equations (21) and (22) are defined in 
terms of cylindrical coordinates 4> and x. In order to apply these equations 
to general wing-body configurations, they must be transformed to local polar 
coordinates. This transformation leads to the following results (see appen- 
dix A for details): 



and 


tan T = cos 6^ 


3$>\ /3*b\ / 3r b\ 

3x/<j>\3$ /x \ 3 X /$ 


(25) 


where 



(26) 




r 0 


Rb sin $ 



(27) 
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( 28 ) 



/3R b \ /3R b \ 

The derivatives and 1 are obtained from the geometry description 

\3* / x \3X /* 

program QUICK. 

Similarly/ equation (19) gives the location of the streamline in terms of 
the cylindrical angle <f> . For wing-body configurations, the angle $ may not 
describe a unique point on the body cross section. (See fig. 4.) To avoid this 
problem, equation (19) can be rewritten in terms of the local polar angle $ 
in the following form (see appendix A for details): 
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3<t> 

and where | — I and 


tively, and 


/3R b \ /3R b \ 

' 3 ® K' W/$' 


are given by equations (19) and (26), respec- 


dr o 

and are calculated in the geometry descrip- 

<3X 


tion program QUICK. 


After obtaining the angle 0, which is related to the velocity direction 
on the surface, from the inviscid solution (appendix B) , equation (30) can be 
integrated to obtain the streamline location on the body surface as a function 
of 4> and x. This is a general equation and can be applied to any configu- 
ration where the pole can be located to give R b as a single- valued function 
of $ and X* (See fig. 4.) For the special case where rg = 0, the local 

polar coordinate system reduces to a cylindrical coordinate system, and equa- 
tion (30) reduces to the form 


/34>\ /3<t>\ 

W/g W/g 


(33) 


as would be expected. 

The metric coefficient h 2 and distance along a streamline s are 
obtained by solving equations (20) and (23) , respectively, as previously 
discussed. 


Calculation of Heating Rates 

Approximate three-dimensional heating rates along surface streamlines can 
be calculated through the axisymmetric analogue using any available axisymmetric 
boundary-layer solution. Although heating rates are usually calculated from a 
solution of the complete axisymmetric boundary-layer equations (eqs. (7) to 
(10)), this is unnecessary because very accurate results can be obtained more 
easily from approximate heat- transfer relations similar to those presented in 
reference 25. Following the nomenclature of reference 25, the expression for 
the heating rate to the surface is given by the following equation (see appen- 
dix C for details) : 


<3w = 


H e (py) w u e h 2 


L N Pr,w 




(34) 
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where the transformed surface distance £ is given by 


-r 


(py) w u e h 2 ds 


( 35 ) 


The parameter is the derivative of the nondimensional enthalpy profile 

and is approximated by the expression (see 


normal to the wall 
appendix C) 


3 (H/H e ) 

. 3n _ 


w 


C w = 0.47< 


/PeUe^ 

1 .48 - 0.23 

\pwPwy 


/p e Pe\ 

\P\MviJ 


1/2 


x (NprC 1 " 0- 24Np r (1 - t e ) 2 ]}^l + 0.11 1 f&J (^aw t W > 

where the velocity gradient parameter 8, the nondimensional adiabatic wall 
enthalpy £ aw , and the recovery factor RF are given by the following . 
expressions: 


(36) 


/du 


e\ 


8 = 


2? 

\ds 


t e (pii) w (u e h 2) 2 


(37) 


? aw = RF + (1 - RF)t e 


(38) 


and 



(39) 


The approximate expression for ^1, (eq, (36) ) has been found to predict values 
that agree with exact similar solutions (tabulated in ref. 26) to within 
±4 percent for the following range of conditions: 
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0 . 1 835 = ^ 1.385 

(40a) 

V^v^w/ 

0.2 = t^ — 1.0 

(40b) 

0.0076 ^ ^ 0.75 

(40c) 

ID 
CO 
VII 
1 cn 
VII 
o 

(40d) 


Boundary-layer edge properties for the heat-transfer calculation are 
obtained by assuming that the flow expands isentropically along an inviscid sur- 
face streamline from the stagnation pressure to the local pressure obtained from 
the inviscid solution. Thus, the properties at the edge of the boundary layer 
(which are assumed equal to the inviscid wall conditions) can be expressed for a 
perfect gas as follows: 


P e = P s (Pe/Ps) 1A 

Pe/ Y 
h ° ■ p-(~ 

and 

u e = \J~2 (H e - h e ) 


(41) 


(42) 


(43) 


The wall temperature is assumed to be known; thus, the nondimens ional 
enthalpy ratio at the wall can be calculated from the equation 


C w = c p T w /H e (44) 

The viscosity ratio V> w A>e is computed from the Sutherland viscosity formula 
(ref. 27) 

y w f 1 + (s/t 6 ) 

HI ” (T w /T e ) + (S/T e ) 
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where S is a dimensional constant defined in the present work as 


S = 110,33 K 


( 46 ) 


RESULTS AND DISCUSSION 

In the sections that follow, the limitations and range of applicability of 
the method will be discussed; also, a comparison of calculated heating rates 
with experimental measurements for a spherically blunted 15° half-angle cone, a 
spherically blunted 80° sweep slab delta wing, and a Space Shuttle Or biter type 
configuration will be presented. 


Limitations and Range of Applicability 

The basic assumption of the axisymmetric analogue is that the cross flow in 
the boundary layer is small (ref. 9) . It has been shown that the cross flow in 
the boundary layer will be small when the streamline curvature on the surface is 
small (ref. 10) or when the wall is highly cooled (ref. 11). The first condi- 
tion is rather restrictive since in general it will be satisfied only for smooth 
bodies (that is, bodies with only small changes in curvature) at small angle of 
attack. The second condition is less restrictive because for many applications 
the wall temperature is only a small fraction of the adiabatic wall temperature. 
However, these conditions are both qualitative and it is necessary to compare 
calculated heating rates with experimental data to validate the theory. 

In the present paper, it is assumed that the Reynolds number is large 
enough for boundary-layer concepts to apply but not large enough to cause tran- 
sition to turbulent flow. The first of these assumptions is necessary since 
the concept of the axisymmetric analogue is based on boundary-layer theory. The 
second assumption is not necessary since Cooke and Hall (ref. 13) have shown 
that the axisymmetric analogue can be extended to turbulent flows. However, 
turbulent flows are not considered in the present paper. 

Further, it has been assumed that the boundary layer follows inviscid sur- 
face streamline paths and that the entropy of the fluid at the edge of the 
boundary layer is constant. This is equivalent to assuming that the boundary 
layer is very thin and that all of the fluid entering the boundary layer has 

passed through the normal portion of the bow shock wave. This is a good assump- 

tion in the nose region and for some distance downstream of the nose for high 
Reynolds number laminar flow. 

In the region far downstream of the nose, the boundary-layer thickness is 
no longer negligible, and the entropy of the fluid entering the boundary layer 
is not constant. In this region, it becomes increasingly important to account 

for boundary-layer growth in calculating the metric coefficient h 2 and vari- 

able entropy in calculating boundary-layer edge conditions, both of which can 
have a strong influence on heating rates. These effects have not been consid- 
ered in the present paper because they would unnecessarily complicate the devel- 
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opment of the method and are considered of small importance for the experimental 
comparisons which follow. These effects can be incorporated in the heating cal- 
culations at some future time. 


Comparisons With Experimental Data 

Spherically blunted 15° half-angle cone .- Calculated heating rates are com- 
pared with experimental measurements from reference 28 on a spherically blunted 
15° half-angle cone at angles of attack of 0°, 10°, and 20°. The experimental 
tests (ref. 28) were performed in air (Y = 1.4) at a free-stream Mach number of 
10.6 and a free-stream unit Reynolds number of 1,31 * 10 6 per meter. The ratio 
of wall to boundary-layer edge stagnation enthalpy C w for these tests was 
0,27. The Prandtl number was assumed to be 0,72, 

The a = 0° case is of particular interest because the flow is axisymmet- 
ric and results obtained by the present theory can be compared with other theo- 
retical results. First, it is well-known that the metric coefficient which 
describes the divergence or convergence of streamlines in axisymmetric boundary- 
layer flows is equal to the body cross-sectional radius (refs. 20 and 29) . To 
check this result, the metric coefficient obtained from the present theory for 
this case is compared with the body cross-sectional radius in figure 9. The 
results for the nose region (x/R n ^ 2) are presented in figure 9(a) and the 
downstream region (x/R n £ 2) in figure 9(b). As would be expected for any axi- 
symmetric body at a = 0°, the calculated metric coefficient h 2 agrees almost 
exactly with the body radius r^. 

The axial distribution of heat- transfer rate for a = 0° is presented in 
figure 10. The present theory is compared with experimental data from refer- 
erence 28 and the results of a nonsimilar axisymmetric boundary layer solution 
computed by the method of reference 21 . Both calculated results are in very 
good agreement with the experimental data. The approximate heat- transfer 
relations used in the present theory (eqs. (34) to (39)) yield results compara- 
ble to those obtained using the differential boundary- layer solution procedure 
(ref. 21 ) . 

The axial heat- transfer distributions along two meridional planes - the 
windward symmetry plane ($ = -90°) and the side meridional plane ($ = 0°) - are 
presented in figure 11 for an angle of attack of 10°, Experimental data from 
reference 28 for two different nose bluntness are presented for comparison. The 
calculated heating rates are in good agreement with the experimental data along 
each meridional plane although there is considerable scatter in the experimental 
measurements in the windward symmetry plane ($ = -90°) for values of x/R n 
greater than about 12. 

The circumferential heating distributions for a = 1 0° are presented in 
figure 1 2 at two axial stations: x/R n = 3.1 and x/R n = 12.0, The abscissa 

for each part of the figure begins in the windward symmetry plane (0 = -90°) 
and continues around to the leeward symmetry plane ($> = 90°) . Heating rates 
calculated by the present method agree well with the experimental data around 
the cone at each axial station. 
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The axial heating-rate distributions at a = 20° are presented in fig- 
ure 13, and the circumferential heating-rate distributions for this angle of 
attack are presented in figure 14. The calculated heating rates are in good 
agreement with the experimental data in each instance. These good comparisons 
are very encouraging because the inviscid surface streamlines for this angle 
of attack (see fig. 15) exhibit large divergence over much of the lower surface 
and side of the cone. This suggests that the small cross-flow analogy can yield 
reasonably accurate heat-transfer results even in regions of relatively large 
streamline divergence for values of C w as large as 0.27. 

In all of the theoretical results presented previously the pole of the 
local polar coordinate system has been coincident' with the body axis. Thus, 
r o (x) is equal to zero; and since the geometry is axisymmetric, SR^/aS is 
also equal to zero. This greatly simplifies the problem since many of the terms 
drop out of the equations and may have the effect of masking errors in the pro- 
gram that would show up later for more complicated geometries. Thus, the invis- 
cid flow-field and heat-transfer calculations were repeated using the shifted 
pole configuration defined by the following expressions: 

For OS x S 1.5, 


r 0 (x) = 0 
For 1 .5 S x S 15.0, 


r 0 (x) = -0.13991 (x - 1 .5) 

which is illustrated in figure 16. The calculations were performed for ot = 0° 
and 10°, and the results are compared with previous heat-transf er calculations 
in figures 17 to 19. The heating rates obtained using the shifted-pole config- 
uration agree closely with those obtained previously when the pole was aligned 
with the body axis (i.e., rQ(x) = 0). With these results, heating rates on 
more complicated body geometries can be computed with added confidence. 

Spherically blunted 80° sweep slab delta wing .- Next, calculated heat- 
transfer coefficients are compared with experimental measurements, from refer- 
ence 30, on a spherically blunted 80° sweep slab delta wing (fig. 20) at angles 
of attack of 0°, 10°, and 20°. The experimental tests in reference 30 were per- 
formed in air (y = 1.4) at a free-stream Mach number of 9.6 and a free-stream 
unit Reynolds number of 3.94 x 10® per meter. The ratio of wall to boundary- 
layer edge stagnation enthalpy C w in the tests was 0.33, and the Prandtl num- 
ber was assumed to be 0.72. 

The axial distribution of heat-transfer coefficients for ot = 0° is pre- 
sented in figure 21 . The heat- transfer coefficients calculated by the present 
theory are seen to be in good agreement with the experimental data. 
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The circumferential distributions of heat- transfer coefficient for a = 0° 
are presented in figure 22 at three stations (l / R n = 4, 8, and 12) measured from 
the beginning of the wing leading edge. These distributions were obtained along 
cross-sectional cuts normal to the wing leading edge (i.e., the dashed lines in 
figure 20) . Even at a = 0° f the heating about a delta wing is not axisymmetric 
because of the geometry; but it is symmetric about the wing leading edge 
($> = 0°) . At the upstream station, 7/R n = 4 ' the present theory underpredicts 
the experimental data slightly around the leading edge; but at the two down- 
stream stations, Z/R n = 8 and the theory is in very good agreement with the 
data. 

The axial distribution of heat- transfer coefficient at a = 1 0° is pre- 
sented in figure 23. As at the lower angle of attack, the heat-transf er coef- 
ficients calculated by the present theory are in very good agreement with the 
experimental data. Circumferential distributions at this angle of attack are 
presented in figure 24, At each station, the heating increases in a direction 
away from the windward symmetry plane ($ = -90°) until it reaches a peak 
slightly ahead of the wing leading edge ($ < 0°) and then decreases rapidly 
around the leading edge on the lee side. The heat-transfer calculations do not 
extend all the way around to the leeward symmetry plane ($ = 90°) because a 
cross- flow shock wave occurs in the inviscid flow field solution on the lee side 
of the body, and the boundary-layer equations cannot be used to accurately com- 
pute the heating in this region. Heat-transfer coefficients predicted by the 
present theory are generally in good agreement with the experimental data. The 
most notable exception occurs near peak heating for 7/R n = however, the 
experimental data point in this region takes an inexplicable dip and may be in 
error . 

The axial distribution of heat- transfer coefficients at a = 20° is pre- 
sented in figure 25. It can be noted that the present theory is in good agree- 
ment with the experimental data. The circumferential distributions at this 
angle of attack are presented in figure 26. The heating at this higher angle of 
attack is qualitatively similar to the heating at a = 10°, except the level on 
the windward side and leading-edge region is much higher. Again, it is found 
that the present theory is in good overall agreement with the experimental data. 

The inviscid surface streamlines for the delta wing at a = 0° are pre- 
sented in figure 27. There is one streamline that lies in the wing leading-edge 
symmetry plane (fig. 27(a)) which divides the streamline pattern into a symmet- 
rical upper and lower part. Near the bottom symmetry plane (fig. 27(b)) where 
the heating is low, the streamlines have very little curvature. However, on the 
wing leading edge (fig. 27(a)) where the heating is relatively high, the stream- 
lines diverge very rapidly. 

The overall agreement between the predicted and measured heating rates on 
the blunt-slab delta wing is very good. This is the first time that three- 
dimensional heating rates on a blunt-slab delta wing have been calculated suc- 
cessfully. This is very encouraging because this configuration is not axisym- 
metric and, although analytic, has many features that are similar to more 
complicated vehicles. 
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Space Shuttle Orbiter type configuration .- Finally, calculated heat- 
transfer coefficients on a Space Shuttle Orbiter type configuration are compared 
with experimental measurements from reference 31 at angles of attack of 20° 
and 25°. The experimental tests in reference 31 were performed in air (y = 1.4) 
at a free-stream Mach number of 7.9 and a free-stream unit Reynolds number of 
1.64 x 1 0 6 per meter. The ratio of wall to boundary- layer edge stagnation 
enthalpy for the tests was 0.31, and the Prandtl number was assumed to be 

0.72. The model used in the experimental tests was 0.0175-scale model of the 
full-scale vehicle shown in figure 28. 

In constructing the mathematical geometry model from the QUICK geometry 
program to use in the inviscid flew field calculations, the canopy was faired 
smooth as shown by the dashed line in figure 28, and the vertical tail and 
reaction jet control pod on the rear of the vehicle were omitted. Difficulties 
were encountered when computing the lee-side flow field with the ’’realistic” 
cross-sectional geometry; thus the lee-side cross section was smoothed out using 
an elliptical segment. (See dashed line in section A-A of fig. 28.) The lower 
side of the vehicle geometry was accurately modeled so that the flow on the 
windward surface could be accurately computed. The modifications to the lee- 
side geometry had no effect on the windward-side flow calculations because in 
all cases the cross-flow velocity went supersonic before the geometry modifica- 
tions were encountered. The inviscid solutions were computed over the forward 
60 percent of the vehicle length, back to where the wing started to flare out 
rapidly. (See dashed line in top view of fig. 28.) At this point, the axial 
velocity at the wing tip became subsonic and the solution could not proceed far- 
ther downstream. (See discussion of limitations of supersonic marching tech- 
niques in refs. 17 and 19.) Thus, the heating rates could only be computed on 
approximately the forward 60 percent of the windward side. 

The axial distribution of heat-transfer coefficient along the windward sym- 
metry plane at a = 20° is presented in figure 29. As can be noted, the heat- 
transfer coefficients calculated by the present theory are in very good agree- 
ment with the ejqserimental data. 

Lateral distributions of heat-transfer coefficient across the lower surface 
of the model at a = 20° are presented in figure 30 for two axial stations: 
x/L =0.4 and x/L = 0.5. The calculated heat-transfer coefficients are in very 
good agreement with the few experimental data points that are available. 

Although no data are available near peak heating, the excellent results noted 
previously for the delta wing in this region lend confidence to the present 
predictions. 

For a = 25°, axial distributions of heat- transfer coefficient are pre- 
sented in figure 31, and lateral distributions are presented in figure 32. As 
at the lower angle of attack, the present theory is in very good agreement with 
the experimental data. 

The inviscid surface streamlines on the windward surface of the Space Shut- 
tle Orbiter type configuration are presented in figure 33. It should be noted 
that near the symmetry plane where the heating rates are relatively low the 
streamlines have little curvature; whereas at the outer edge of the body where 
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the heating rates increase sharply, the streamline curvature also increases 
rapidly - as would be expected. 

Thus, on the windward surface of a Space Shuttle Orbiter type configuration 
where the inviscid flow field can be accurately computed the present theory has 
been shown to predict heating rates that are in very good overall agreement with 
experimental data at moderately high angles of attack. 

From the present study, it appears that accurate three-dimensional heating 
rates can be computed on most regions of a vehicle where boundary-layer concepts 
apply if £ w is not too large (i.e., less than approximately 0.4) and an accu- 
rate inviscid flow field solution can be computed. The range of £ w might be 
extended even higher through future comparisons with other experimental data. 
From the experience gained in the present study it is estimated that the heating 
over a "complete" vehicle can be obtained in a few minutes computing time on a 
Control Data CYBER 1 75 (or equivalent) computer once the inviscid flow field has 
been obtained. Consequently, the present theory should prove very useful in 
studying the heating on advanced Earth entry vehicles. 


CONCLUDING REMARKS 

A theoretical method has been developed for computing approximate laminar 
heating rates on three-dimensional conf igurations at angle of attack. The 
method is based on the axisymmetric analogue which is used to reduce the 
three-dimensional boundary-layer equations along streamlines to an equivalent 
axisymmetric form by using the metric coefficient which describes streamline 
divergence (or convergence). The method has been coupled with a three- 
dimensional inviscid flow field program for computing surface streamline paths, 
metric coefficients, and boundary- layer edge conditions. Using this method, 
accurate laminar heating rates can be computed for a wide range of three- 
dimensional configurations at moderately large angles of attack. This conclu- 
sion is supported by good comparisons with experimental data on a spherically 
blunted 15° half-angle cone, a spherically blunted 80° sweep slab delta wing, 
and a Space Shuttle Orbiter type configuration at angles of attack up to 25°. 
Computations proceed very rapidly with the heating calculations for a complete 
configuration requiring only a few minutes of computing time. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
July 1 5, 1 980 
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APPENDIX A 


TRANSFORMATION FROM CYLINDRICAL COORDINATES 
TO LOCAL POLAR COORDINATES 

The equations relating local polar coordinates and cylindrical coordinates 
are given by equations (17a) to (17d) 


r b = Rb\/^ + l cos $> 


<A1) 


4> = tan -1 (A) 


(A2) 


x = x 


(A3) 


where 


A = tan $ + rg(x)/(Rb cos $) 


(A4) 


The transformation operators for transforming derivatives from cylindrical 
coordinates on a surface (4>,x) to local polar coordinates on a surface ($,x) 
are 



/3$\ /3 \ /3y\ /a \ 

\34>/ X \W x \ 9< ty X \ 9 X/ $ 


(A5) 


i 3 \ /3$\ /3 \ /3x\ /3 \ 

V3x/<J) WcfrWx 


Now since 



and 
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equations (A5) and (A6) can be written as 



(A7) 


and 


/3 \ / 3 ®\ /: 3 \ /3 \ 

\3x/ (jj \3x/$ \<My x \%X/$ 


Similarly, the reverse transformation operators are 



(A9) 


3 \ /a \ /a \ 

3X/$ \9xA Wx \3x/<J) 


(A10) 


Now solving equation (A9) for 



the following result is obtained: 



(All) 


Comparison of equations (A7) and (All) shows that 
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(A12) 


Similarly, solving equation (A10) for 


replace 



results in 



and using equation (All) to 



(A13) 


Comparing equation (A10) with equation (AT 3) gives 



(AT 4) 




/3<fr\ 

The derivatives 

— and 

1 ) 

follows : 

Wx 

\ 9 X/$ 


can be obtained from equation (A2) as 



(A15) 


or 
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r 3<t>\ _ 1 /3 a\ 

3$ /X 1 + A 2 W/X 


(AT 6) 


Similarly, 


3<t>\ _ 1 /3 a\ 

3 x/$> i + a 2 W$ 


(A17) 


From equation (A4) , the following expressions can be obtained for the deriva- 
( 3a\ /3a\ 

tives — I and — , respectively; 

\3*/ x \3X/« 


'3A 


^Rb 


, R 2 + rnfRh sin ® - 1 1 cos $ 

34>/v t >2 cos 2 °L b 


X R£ 


(AT 8) 


/3a\ 

1 


dro 

r o/ 3R b\ ~ 

\ax/« 

R b cos 

$ 

dX 

Rb W / $ 


(AT 9) 


/3<t> \ 

Combining equations (AT 6) and (A18) gives the following expression for — : 

W 


raj' 

^3#> 


R b + r 


f (**b\ 

0 (R b sin $ - y 


cos $ 


(l + A 2 ) (R b oos 2 $) 


(A20) 


Similarly, from equations (A17) and (AT 9), 


27 



APPENDIX A 


dr 0 r 0 /9Rb'| 

'34>\ W ~ Rb\5x/ $ 

^ + j2'j cos $) 


From equations (A12) and (A20) , 



Similarly/ from equations (A14), (A20), and (A21 ) , 



-R b cos 


r 


dr 0 r 0 /3R b \ 

dX R bW ) <I> 


| 2 f /3R b \ 

*b + r 0 |R b sin * - j 


COS $ 


(A21) 


(A22) 


(A23) 


In summary, the transformation operators defined by equations (A7) and (A8) can 

/3$\ 

be used with equation (A22) which defines 

'a$\ 




and equation (A23) which 


defines — j to transform derivatives from cylindrical coordinates to local 


^ x /4> 


dr. 


'3 R b ' 


'3Rh 


polar coordinates. The quantities Rb, rg, , ( I , and I 


dX 


l 3$ 


^X h 


used 


in the previous equations must be obtained from the geometry description 
program. 


These operators will now be used to transform the equations for the body 
geometry angles 6^ and T from cylindrical coordinates to local polar 
coordinates. First, consider 6^ which is given by equation (21) 
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1 / 3r b\ 
tan 5* = — — / 
*bW / x 


Using equation (A7) , the derivative j can be written as 

\ 3< t>/x 


3r b\ _ /a#\ / 3r b\ 

9T/x = WxWx 


Fran equation (A1 ) , 


3r b\ | lo 1 / 3R b\ A /8 a\ 

= RbVT + A cos $ — + — - tan $ 

8 »/ x ^ *»V»/ X , ( jiUi, 


Substituting this result into equation (A24) and using equation (AT) gives 


tan 6* = 


3®\ / 3r b\ 


R bv + °° s ® 


where 

/a# 

(A26) . 



Next, consider T which is given by equation (22) 

/ 3r b\ 

tan T = cos 6* 


(A28) 
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(**b\ 

From equation (A8), the derivative 1 can be written as 

\3*/<i> 


/ 9r b\ /a*\ / 9r b\ / 9r b\ 

W/<j> ~ W<AW X + \ix/# 

From equation (AT ) , 



= Rb'A 


+ A 2 


cos $ 


l_/9Rb\ 
R b\ 9 X U 


+ 



1 + A 2 '^/* 


(A29) 


(A30) 


Substituting equation (A29) into equation (A28) yields 


tan T = cos <S 


<J> 


3$\ / 9r b\ / 9r b\ 

9x / 4> \ 9 ® /x / $ 


(A31 ) 


/3*\ / 9r b\ 

where — is given by equation (A23) , I 1 is given by equation (A26) , 


V 9x /<t» 

/ 9r b\ 

and is given by equation (A29) . 

\ 9 X / $ 


v3$ 


Next, consider the transformation of the equation describing the streamline 
location (eq. (19)): 



1 /tan 0 cos 6^ 
r^V cos T 


- tan T sin 



(A3 2) 


from cylindrical coordinates 4>,x to local polar coordinates $,X» Using 
equations (A2) and (A4), the following are obtained: 


'3<t>\ 1 /3 a\ 

? x ) & ! + a 2 We 


(A33) 
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and 


'3a\ 

W3 


I R b + r 0 sin 1 

R b cos 2 $ / \3 x/ 3 R b cos <J> 


'3r 0 \ r 0 /3R b \ 

^3x ) $ R b \3x/3 


(A34) 


/3$\ 

Now, combining equations (A33) and (A34) and solving for — yields 

W3 


^3$\ /3$ 


^X/B \ 9x /B 


SL - (’ ♦ * 2 )l 


Rfc> oos 2 9 \jfo\ 
R b + r 0 sin « A\3x/3 


Rfc COS $ 


r o/ 9R b\ / 9r °'\ 
R b\ 9x /3 \ 9x /0 


(A3 5) 


Next, changing from an x, 3 coordinate system to an x,4> coordinate system. 


I 9 I \ \ / 9 < t>\ f 9 \ 

l 9x /3 \ 9 x / 4 > \3x/3\3<t>/ x 


(A36) 


Further, using equations (A7) and (A8) , 


- 

V 9x / 3 


V 9x /3 


'3®\ /3 \ 

, 9 <VxWx 


^3$ \ /3 \ /a \ 

V 9x / <j>\3 v X \ 9 X/ $ 


(A3 7) 


or rearranging terms. 


Us 


3 <}) \ / 34>\ /30 

9x /3\ 9< 1 Vx \ 9x /<j> 


f 9 A / 9 

lii/v + ~~ 




(A38) 


Now, since ro = ro(X)» 
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f3r 0 \ /3r 0 \ o 

^ x /e \3X/# dx 


(A39) 


With the use of equation (A38) , the following is obtained: 


/'3Rb\ 

<3x /g 


/3<j> \ /3«\ 

v3 x /<f> \3x/g \34>/x 


W\ /3Rb\ 
^3$ / x \3x/$ 


(A40) 


/a*\ M 

where — I is given by equation (A23), — I is given by equation (A32) , 


W 


4> 


We 


3*\ / 3R b\ f dR b) 

is given by equation (A22) , and I 1 and I ) are given by 


the geometry description program QUICK. 


.3$ 


k 3 X/ $ 
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CALCULATION OF SURFACE VELOCITY DIRECTION FROM 
INVISCID SOLUTION 

In order to integrate equation (30) to determine the location of a stream- 
line on the surface, it is necessary , to know 0= 0( x) • This must be obtained 
from a solution of the inviscid flow field. The solution is obtained using the 
STEIN (supersonic three-dimensional external inviscid) flow field code described 
in references 17 and 19, From the inviscid solution the Cartesian velocity com- 
ponents are known; these can then be used to determine the velocity direction 0 
on the surface. 

Consider sketch (a), which shows the Cartesian velocity components w*| , W 2 , 
and W 3 , and sketch (b) , which shows a cross-section plane normal to the body 
axis: 



Sketch (a) Sketch (b) 


The velocity components v t and v n are given by the equations 

v t = W3 cos (4) - S<j,) - w 2 sin (<(> - 6 ^) (B1 ) 

v n = »3 sin (4> ~ 6 <t)) + w 2 cos (<(> - 6 ^) (B2) 
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Now, in a plane parallel to the axis and normal to the body shown in sketch (c). 


v 



Sketch (c) 

the velocity component tangent to the body w t is given by the equation 
w t = v n s ^ n T + w-j cos r 

Finally, the streamline direction 0 is given by the equation 
0 = tan - ^ (v t /w t ) 


(B3) 


(B4) 
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DEVELOPMENT OF EQUATIONS FOR CALCULATING 
HEAT-TRANSFER RATE AT WALL 

The heat-transfer rate to the wall q w is given by the expression 


<3w = 



(Cl) 


3T 

where k is the thermal conductivity of the fluid and I — I is the derivative 

\3n/ w 

of the temperature of the fluid normal to the wall. Equation (Cl) can be 
rewritten in terms of the static enthalpy derivative as follows: 


k w /3h\ 
C P, w vW w 


(C2) 


Neglecting cross flow in the boundary layer, the static enthalpy h is related 
to the total enthalpy H by the equation 


h = H 


u 2 + w 2 
2 


(C3) 


Differentiating the terms in equation (C3) with respect to n yields 


3h 3 h 3u 3w 

3n 3n 3n 3n 


and since u = w = 0 at the wall. 


(-) - ( 2 ) 

\3n/ w V3n/ W 


(C4) 
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Substituting equation (C4) into (C2) yields 


<3w = 


k w 


c P»w\ 3n /w 


(C5) 


Multiplying and dividing the right-hand side of equation (C5) by y w yields 


~ ^w 


f k \ /3 h\ ^w /3h\ 
\. c pVj w\ 3n / w N Pr, w\ 3n / w 


(C6) 


The laminar boundary- layer equations are usually solved in a Levy-Lee type 
coordinate system (£,f)) given by the following equations (ref. 25): 


-r 


(PU ) W u e h 2 ds 


(C7) 


u e h2 pn 

T] = IT X 


p dn 


(C8) 


Using these variables, equation (C6) can be transformed to the following form: 


q „ . (Z 


(py) w u e h 2 /9 H \ 


N Pr, w\ 3n ' w\ 3 ^/w \^Vw 

N Pr, w» 2 5 


(C9) 


The total enthalpy in equation (C9) can be nondimens ionali zed by multiplying and 
dividing the right-hand side by H e (which is independent of n ) to obtain the 
following: 


q w 


H e (PU) w u e h 2 


L N Pr 


.M 


3 (H/H e ) 

. 3n _ 


»w 


(CIO) 
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Now defining 


C 


i 

w 


3 (H/H e ) 


L 


■J w 


(Cll) 


equation (CIO) can be written as 


q w = 


H 0 ( Py ) w u e h 2 


L Np r , w V2 


»w 


(Cl 2) 


which is the same as equation (34) . 

In reference 25 it was shown that could be approximated by the 

expression 


0.475 


= 0.47 


/ p eM e y . . 

— [i - n - N PrrW )(i - t e )] 

\PwM w / 


x (1 + o.i\/Bt^)(i - c w ) 


where 3 is the velocity gradient parameter given by the equation 


(Cl 3) 



t e (PU) w (u e h 2 ) 2 


(Cl 4) 


Equation (Cl 3) predicts values of that agree with tabulated results from 

boundary- layer solutions presented in reference 26 to within ±1 0 percent for 
the following range of conditions: 


/ Pe p e\ 

0.1835 S £ 0.9367 

\p w y w / 
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0.2 -< t e ^ 1.0 
0.0076 ^ C w ^ 0.75 
0.5=6= 3.5 


Thus, from equations (Cl 2) and (Cl 3), heating rates can be rapidly calculated 
from boundary-layer edge and wall parameters without having to solve the com- 
plete set of boundary-layer partial differential equations (7) to (10). 

In the present study, it was found that a slightly more accurate prediction 
of could be obtained by using the following equation (same as eq. (36) ) : 



where C aw is given by 

S aw = RF + (1 - RF)t e (Cl 6) 

and the wall recovery factor RF for laminar flow is assumed to be 


*2 “ \/ N Pr,w 


(Cl 7) 


Predicted values of Cw from equation (Cl 5) agree with results of the boundary- 
layer solutions presented in reference 26 to within approximately ±4 percent for 
the following range of conditions: 


/PeM 

0.1835 S ^ 1.385 (Cl 8a) 

\P w^w/ 


0.2 £ t e ^ 1.0 (Cl 8b) 

0.0076 ^ ^ S 0.75 (Cl 8c) 
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0 £ 3 S 3.5 


(Cl 8d) 


A few typical comparisons are presented in table I. 

The range of conditions over which equation (Cl 5) has been correlated 
covers the range of experimental conditions of interest in the present study 
which are given for experimental data by the following: 


/Pe^eN 

0.6 £ ^1.1 (Cl 9a) 

\P w Pw/ 

0.4 ^ t e ^ 1.0 (Cl 9b) 

0.27 S £ w ^ 0.33 (Cl 9c) 

0 ^ 3 ^ 0.7 (Cl 9d) 


Equations (Cl 2) (or eq. (34)) and (Cl 5) (or eq. (36)) have been used in the 
present study to provide a simple but accurate method of calculating heating 
rates along a surface streamline. 
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TABLE I.- COMPARISON OF VALUES OF C w FROM EQUATION (Cl 5) WITH EXACT 
SIMILAR BOUNDARY-LAYER SOLUTIONS FROM REFERENCE 26 



B 

£ w 

Npr 

l 


Percent 

error 

PwMw 

Reference 26 

Equation (Cl 5) 

0. 51 22 

1.0 

0.1 52 

0.709 

0 

0.2845 

0.2755 

-3.1 6 

1 

1 

1 


.5 

.3055 

.2969 

-2.82 

T 

\ 

\ 

t 

3.5 

.3377 

.3322 

-1 .63 

.9149 


.75 

.735 

0 

.101 3 

.1049 

3.55 

1 


1 

1 

.5 

.11 41 

.1130 

-.96 

T 

\ 

\ 

T 

3.5 

.131 6 

.1265 

-3.88 

1.3850 

.2 

.5 

.699 

0 

.1 518 

.1479 

-2.57 

i 


1 

1 

1.0 

.1697 

.1641 

-3.30 

T 

1 

T 


3.4 

.1 817 

.1 778 

-2.1 5 

. 2322 

no 

.01 52 

.709 

0 

.2230 

.2256 

1.17 

| 

I 

1 

1 

.5 

.2399 

.2431 

1.33 

T 

\ 

T 

f 

3.4 

. 2659 

.271 3 

2.03 

.2505 

.8 

. 01 52 

.709 

0 

.2206 

.2246 

1 .81 

1 

| 

| 

1 

.5 

.2365 

.2421 

2.37 

\ 

\ 

T 

Y 

3.2 

.2616 

.2688 

2.75 

. 3462 

.6 

.03 

.680 

0 

.2293 

.2308 

.65 

1 

1 

1 

1 

.5 

.2454 

.2487 

1.34 

T 

\ 

\ 

Y 

3.0 

.271 5 

. 2747 

1.18 

.6249 

.4 

.10 

.768 

0 

.2891 

.2940 

1 .69 

1 

1 

1 

1 

.5 

. 3099 

.31 69 

2. 26 

T 

\ 

\ 

1 

3.5 

.3498 

.3545 

1.34 

. 3897 

.8 

.05 

.768 

0 

.2871 

.2877 

.21 

1 

1 

1 

1 

.5 

.3087 

. 31 01 

.45 

T 

\ 

| 

\ 

3.2 

. 3426 

.3444 

.53 

.7931 

.2 

.10 

.768 

0 

.2889 

.2947 

2. 01 

1 

1 

1 

I 

.5 

. 3086 

.31 77 

2.95 

\ 

1 

\ 

\ 

3.5 

.3518 

.3554 

1.02 

.9367 

.6 

.50 

.699 

0 

.1723 

.1700 

-1.33 

1 

1 

1 

1 

.5 

.1882 

.1 832 

-2.66 

| 
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Y 

3.5 

.2118 

.2050 

-3.21 

1.0810 
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. 50 

.699 
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.1 628 
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1 
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.1763 
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1 

\ 
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Figure 6.- Streamline pattern of nose of general three-dimensional vehicle. 











Figure 10.- Axial distribution of heating 
cone at a = 0°. = 10.6; Y = 1*4; 

Rn = 0.95 cm; q s = 224.32 kW/m^-s. 



:ate on a spherically blunted 15° half-angle 
N Re f oo = 1*31 x 10® per meter; £ w = 0.27; 



Figure 11.- Axial distribution of heating rate on a spherically blunted 15° half-angle 
cone at a = 10°. = 10.6; y = 1 .4; N R6f0o =1.31 * 10 6 per meter; £ w = 0.27; 

q s = 224.32 kW/m 2 -s. 
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(ci) x/R n ~ 3.1 . (b) x/R n - 12.0. 

Figure 12.- Circumferential distribution of heating rate on a spherically blunted 15° 
half-angle cone at a = 10°. ■ 10.6; Y = 1.4; NR e/00 = 1 .31 x 10° per meter; 

= 0.27; q s = 224.32 kw/m 2 -s. 
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Figure 13.- Axial distribution of heating rate on a spherically 
blunted 15° half-angle cone at a = 20°. * 10.6; Y = 1 .4 

%e oo = 1 * 31 x 106 Pe r meter; £ w » 0.27; q s ■ 224.32 kw/m 2 
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(a) x/R n = 3.1 . 


(b) x/R n = 12.0, 


Figure 14.- Circumferential distribution of heating rate on a spherically blunted 15° half 
angle cone at a = 20°. M = 10.6; Y = 1.4; N R6f00 = 1 .31 x 1 0 6 per meter; S w = 0.2 
q s = 224.32 kW/m 2 -s. W 



5 Side view of inviscid surface streamline P^ tern 
5 " re J';T ® iv blunted 150 half-angle cone at a - 20 . 
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Figure 16.- Geometry of shifted pole for spherically blunted 15° half-angle cone. 



Figure 17.- Comparison of axial heat-transfer calculations on a spherically blunted 15° 
half-angle cone with shifted and unshifted poles at a = 0°. = 10.6; Y = 1 •*; 

N Re,<*, x 106 Pe r meter; C w = 0.27; R n = 0.95 cm; q s = 224.32 kW/m 2 -s. 
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Figure 18.- Comparison of axial heat-transfer calculations on a spherically blunted 15° 
half-angle cone with shifted and unshifted poles at a = 10°. = 10.6; y = 1 .4; 

N Re,oo = 1 *31 x 10® per meter; C w * 0.27; R n - 0.95 cm; q s = 224.32 kW/m^-s. 


(a) x/R n =3.1. 


(b) x/R,, = 12.0. 


Figure 19.- Comparison of circumferential heat-transfer calculations on a spherically blunted 15° 
half-angle cone with shifted and unshifted poles at ot = 10°. = 10.6; Y =1.4; 

Nrp m ■ 1.31 x 10 6 per meter; C w = 0.27; R n = 0.95 cm; q,, = 224.32 kW/m 2 -s. 




Figure 20.- 80° sweep slab delta wing. R n = 1.27 cm. 






(a) 2/R n = 4. (b) 2/R n = 8. (c) 2/R n = 12. 

Figure 22.- Circumferential distribution of heat-transfer coefficient on a spherically blunted 
80° sweep delta wing at a - 0°. = 9.6; y = 1.4; N Re = 3.94 x 10 6 per meter; 

C w = 0.33; R n = 1.27 cm; ho = 0.3796 kW/m 2 -s-K. 
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Figure 23.- Axial distribution of heat-transfer coefficient in windward symmetry 
plane of spherically blunted 80° sweep slab delta wing at a = 10°. = 9.6 

Y = 1.4; N R6f0O = 3.94 x 10® per meter; C w = 0.33; R n = 1.27 cm; 
h q = 0.3796 kW/m 2 -s-k. 
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Figure 25.- Axial distribution of heat-transfer coefficient in windward 
symmetry plane of a spherically blunted 80° sweep delta wing at 
a = 20°. ^ = 9.6; y = 1 .4; NRe r oo = 3.94 x 10® per meter; 

= 0.33; R n = 1.27 cm; h c = 0.3796 kW/m 2 -s-K. 
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Figure 26.- Circumferential distribution of heat-transfer coefficient on a spherically blunted 
80° sweep slab delta wing at a = 20°. = 9.6; y = 1.4; N R6f00 = 3.94 x 10 6 per meter; 

5 W = 0.33; R n = 1.27 cm; h s = 0.3796 kW/m 2 -s-K. 
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x/L =0.4. 


x/L = 0.! 


Figure 30.- Lateral distribution of heat-transfer coefficient 
on Space Shuttle Orbiter type configuration at a = 20°. 

M ro = 7.9; Y = 1 .4; % e ,oo = 1 * 64 x 106 per meter; 

C w = 0.31; h r = 0.4984 kW/m^-s-K. 
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Figure 31 Axial distribution of heat-transfer coefficient 
on Space Shuttle Orbiter type configuration at a = 25°. 
M oo = Y = 1*4; N Re,oo = 1 *64 x 10® per meter; 

C w = 0.31; h r = 0.4984 kW/m 2 -s-K. 


1. Report No. 

2. Government Accession No. 

NASA TP-1698 


4. Title and Subtitle 


CALCULATION OP LAMINAR HEATING RATES ON 

THREE-DIMENSIONAL CONFIGURATIONS USING THE 

AXISYMMETRIC ANALOGUE 


7. Author(s) 


H. Harris Hamilton II 


9. Performing Organization Name and Address 


NASA Langley Research Center 

Hampton, VA 23665 


12. Sponsoring Agency Name and Address 


National Aeronautics and Space Administration [ 

Washington, DC 20546 



3. Recipient's Catalog No. 


5. Report Date 

September 1980 


6. Performing Organization Code 


8. Performing Organization Report No. 

L-13748 


10. Work Unit No. 

506-51 -1 3*07 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Paper 


14. Sponsoring Agency Code 


15. Supplementary Notes 

The information presented in this report was included in a report submitted in 
partial fulfillment of the requirements for the Degree of Engineer, The George 
Washington University, December 1979. 


16. Abstract 


A theoretical method has been developed for computing approximate laminar heating 
rates on three-dimensional configurations at angle of attack. The method is 
based on the axisymmetric analogue which is used to reduce the three-dimensional 
boundary- layer equations along surface streamlines to an equivalent axisymmetric 
form by using the metric coefficient which describes streamline divergence (or 
convergence) . The method has been coupled with a three-dimensional in viscid flow 
field program for computing surface streamline paths, metric coefficients, and 
boundary- layer edge conditions.. 


17. Key Words (Suggested by Author(s)) 

Axisymmetric analogue 
Three-dimensional heating 
Boundary- layer flow 


18. Distribution Statement 


Unclassified - Unlimited 


19. Security Classif. (of this report) 
Unci ass if ied 


20. Security Classif. (of this page) 
Unclassif ied 


Subject Cateqor 


21. No. of Pages 22. Price 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1980 





















National Aeronautics and 
Space Administration 

Washington, D.C. 

20546 

Official Business 

Penalty for Private Use, $300 



10 1 10, D, 090380 S009Q3DS 

DEPT OF THE A IE FORCE 

IF NEAPONS LABOBATOBX 

ATTN; TECHNICAi LIBBAEX (SUL) 

KIHTLAND AFB NM 87117 


Postage and Fees Paidi*-^ f 
National Aeronautka|iOd^] 
Space Administratiom^^ 
NASA-451 



I, V:. 


NASA 





